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Abstract Marginal likelihoods for the cosmic expan- 
sion rates are evaluated using the 'Constitution' data 
of 397 supernovas, thereby updating the results in some 
previous works. Even when beginning with a very 
strong prior probability that favors an accelerated ex- 
pansion, we obtain a marginal likelihood for the decel- 
eration parameter qo peaked around zero in the spa- 
tially flat case. It is also found that the new data sig- 
nificantly constrains the cosmographic expansion rates, 
when compared to the previous analyses. These results 
may strongly depend on the Gaussian prior probabil- 
ity distribution chosen for the Hubble parameter rep- 
resented by h, with h — 0.68 ± 0.06. This and simi- 
lar priors for other expansion rates were deduced from 
previous data. Here again we perform the Bayesian 
model-independent analysis in which the scale factor is 
expanded into a Taylor series in time about the present 
epoch. Unlike such Taylor expansions in terms of red- 
\ shift, this approach has no convergence problem. 
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1 Introduction 

It is generally accepted that a more appropriate way to 
measure the acceleration of expansion of the universe 
is to resort to a cosmographic or model-independent 
analysis. In the conventional model-based analyses of 
distance modulus-redshift {jJi—z) data of Type la super- 
nova (SN la), the accelerated expansion of the universe 
is an indirect inference based on the best fit values of pa- 
rameters, such as the density parameters Q m , etc. 



On the other hand, in a model-independent approach, 
the scale factor a(t) is expanded as a Taylor series in 
time about the present epoch and the marginal likeli- 
hoods of its coefficients are computed using the data. 
The marginal likelihood for the deceleration parameter 
gives an estimate of the acceleration of cosmic expan- 
sion. Since practically one has to truncate the series to 
some finite order, the basic assumption here is that a(t) 
is expressible as a truncated Taylor series or polyno- 
mial. Evaluating the deceleration parameter by adopt- 
ing this method, it was confirmed model-independently 
that the universe is un dergoing an accelerated expan- 
sion |johnll2004il2005l ). 



In this paper we report the updating of the marginal 
likelihood for each of the expansion coefficients found 
in the above work. This is performed for the case 
of a fifth order polynomial. A notable result in the 
present Bayesian model-independent analysis is that 
even when beginning with a very strong prior probabil- 
ity that favors an accelerated expansion, the marginal 
likelihood for the deceleration parameter qo is found 
peaked around go = in the spatially flat case. It is 
also found that the new data significantly constrains the 
cosmic expansion rates appearing in the Taylor expan- 
sion, when compared to the previous data. We also note 
that successive terms in the series decreases sufficiently 
fast, thereby verifying the assumption of a converging 
Taylor series in time for the cosmic scale factor. 

Other model- independent approaches, which Taylor 
expand the distance modulus /i in terms of redshift 



eg. ( Shapiro & Turner 


2006; Cattoen & Visser 


2007 


Guimaraes et al. 20091; 


Seikel & Schwarz 2009)] 
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Guimaraes et al. II2009I ). The argument behind this as- 
sertion is as follows: For an expanding universe, z < 
corresponds to the future and z = — 1 is the redshift 
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when the universe has expanded to infinite size. Since 
z = — 1 is a pole, by standard complex variable theory, 
the radius of convergence of a series about z = is 
atmost | z |= 1, so that it fails to converge for z > 1. 
When compared to this, our approach of expanding the 
scale factor in terms of t about the present epoch to 
is advantageous, for the series converges for all times. 
Even the lookback time T = t — to is evaluated by nu- 
merically solving an equation which involves a Taylor 
scries in time. Hence there is no convergence problem 
in the present work. However, it should be noted that 
all analyses which make use of such Taylor expansions, 
in practice, employ polynomials and hence convergence 
is not a serious problem. For instance, one can see that 
there is convergence in certain special cases of t he low 



order polynomial fit by iGuimaraes et al. I ((2009) 



2 Marginal likelihoods for the cosmic 
expansion rates 

With t — t = T, where to is the present time, the scale 
factor of the universe is expan ded into a Taylor series 
about the present epoch to as ( John 20041 2005 ) 



used to find T\ in terms of z, for each combination of 
parameter values. This is done in a direct and purely 
numerical way. 

We may now obtain the distance modulus as fi = 
5 log (D/lMpc) + 25. Here D and hence \x are functions 
of z and contain parameters h, qo, ro, so, uo, etc. In 
the following, we keep only terms up to fifth order in 
the Taylor expansion. 

The likelihood function is 

C = cxp[~x 2 (h, qo,r , so, u Q )/2] 
where x 2 is given by 



2 { £fc~ (J>k(zk;h,qo,r :) 
X = 



(5) 



Here (Ik is the measured value of the distance modu- 
lus of the k th supernova, fik(zk', h, qo, • ■) is its expected 
value (from theory) and Ofc is the uncertainty in the 
measurement. 

The likelihood for the truncated Taylor series form 
of scale factor can be found as 



a(t + T) = a x 
1 + H T - + - ££My4 + ^i T s 4 

Note that T assumes negative values for the past. For 
a k = flat universe, one can put ao = 1. The remain- 
ing parameters in the theory are the Hubble parameter 
Ho = lOO/i km s" 1 Mpc— 1 , the deceleration param- 
eter qo and higher order expansion rates such as ro, 
s , u , etc. Our task is to deduce the values of these 
parameters from the observational data. 

For a light pulse emitted from a SN situated at the 
coordinate r\ at time t\ and reaching us at r = at 
time to, the RW metric allows one to write 



(l)C(M) = J dhj dq J dr J ds J du 

p{h)p(q )p(r )p(so)p(u ) e" x2/2 . (6) 

where p(h)p(qo)p(ro)p(so)p(uo) is a product of Gaus- 
sian probability distributions of each of the parameters. 

This is an approximation to p(h,qo,ro ), the prior 

probability distribution to be introduced in equation 
©■ 

The marginal likelihood of any one parameter can be 
computed by integrating the likelihood function, multi- 
plied by an appropriate prior probability distribution, 
over all parameters except the concerned one. For in- 
stance, the marginal likelihood for go can be found as 



fit, 



cdt 
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dr 



, tl u W Jri (l-fcr 2 )V2- 
For a k = RW metric, this can be used to obtain 



(2) 
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(3) 



With this, we can compute the luminosity distance D = 
riao(l + z). An important part of the calculation is the 
solution of the equation 



l + z = 



a(t ) 



a(to + Ti)' 



(4) 



£(<7o) = / dh J dr J ds J du 
p(h)p(r )p(s )p(u ) e~ x2/2 . 



(7) 



A salient feature in the present computation is that 
while using ([7J, the marginal like l ihood s obtained in 
the previous analysis ( John 2004L 2005) are taken as 
the prior probability distributions, for the correspond- 
ing coefficients. These references have used flat priors, 
since there were no other previous work evaluating these 
marginal likelihoods. But it was proposed there itself 
that the posterior marginal likelihoods obtained in it 
shall be used as priors for future analyses. We note 
that the present work is the appropriate place to make 
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use of this. However, it would not be computationally 
feasible to use the posterior in the previous analysis as 
prior, in terms of a table of values. Therefore, we ap- 
proximate those distributions by Gaussian functions, 
with the cor responding m ean and standard deviations 



obtained in (jJohn 1 120051 ) . A comparison with the ac- 



tual plots show that this is a reasonable approximation 
for most coefficients. The product of such individual 
priors is the combined prior, used in equation (J7J. 

It should be noted that the ranges of flat priors in the 
previous analysis wer e chosen on the basis of the same 
'All SCP' SN data in lKnop et. al~l (|2003l) itself, even- 



though it ran the risk of using the same data twice. 
There we first found the ranges of the 'contributing' 
values of the parameters by varying them arbitrarily. 
Later, flat prior probabilities were assigned for those 
ranges and they were used to find the posterior marginal 
likelihoods. In the present analysis, the marginal like- 
lihoods thus obtained are used as priors. 

As mentioned above, the present mo del-independent 
analy sis uses the 'Constitution' data (jHicken et. al 
2009) of 397 SN. In this connection, it shall be noted 



that some SN la, which appeared in the dataset of 54 
'All SCP' SN used in the previous analysis are present 
in the Constitution set too. But the values of fik, Zk 
and the errorbars of such objects are found modified 
to some extent in the new release. Therefore, it was 
opted not to exclude such SN from the Constitution 
set. We note that at any rate, the Gaussian priors ob- 
tained from the previous analysis are a better option 
than flat priors. 

We have computed the marginal likelihoods of four 
important expansion rates of the present universe, 
namely qq, ro, sq and uo, and the results are shown 
in Figs. (l)-(4). Terms up to fifth order are kept in the 
expansion, but only the flat (k — 0) case is considered. 
This is equivalent to assuming a ^-function prior for 
the flat spatial geometry. The joint prior probability 
used for other parameters was the product of individ- 
ual Gaussian functions in each parameter with mean 
and standard deviations as follows: h = 0.68 ± 0.06, 
q = -0.90 ± 0.65, m =2.7 ±6 .7, s =36.5 ±52.9, and 
u Q 



142.7 ± 320 (|johnll2005l) . In each case, the inte- 



grations were performed in the 2a range of each of the 
parameters. We have performed variation with respect 
to h, though marginal likelihood for this parameter was 
not drawn. The step sizes chosen for these parameters 
were Ah = 0.01, Aq a = 0.1, Ar = 1, As = 20 and 
Au = 100. 

The results show that there is significant constrain- 
ing of the parameters while using the new and refined 



data, compared to the corresponding results in (jJohn 



Fig. 1 Marginal likelihood for the parameter go (in units 
of 10~ ), while using the polynomial of order 5 




2005). It is to be reminded that the marginal likeli- 
hoods are not precisely probability distributions for the 



Fig. 2 Marginal likelihood for the parameter ro (in units 
of KT 105 ), while using the polynomial of order 5 



parameters; instead, they are the probability for the 
data, given the model and the parameter values. How- 
ever, we here compute mean and standard deviations, 
considering the marginal likelihoods as distributions. 
The new mean and standard deviations are the follow- 
ing: q = 0.04±0.30, r = -4.5±4.6, s Q = -42.8±52.5, 
and uq = 320. 5 ± 2 1 3.0. T he marginal likelihood for qo 
obtained by I John I (|2005t ) , which is also used as prior 
for this parameter in the present work, is reproduced 
here in Fig. ([5]) for comparison. The standard devi- 
ations of each of these parameters, except that of sq, 
have decreased substantially and this leads to the above 
assertion that the Constitution data constrains the cos- 
mic expansion rates significantly. 

It shall be noted that even when beginning with 
a prior probability distribution centred around qo = 
—0.9, which is strongly in favor of an accelerated ex- 
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Fig. 3 Marginal likelihood for the parameter so (in units 
of 1(T 105 ), while using the polynomial of order 5 





Fig. 4 Marginal likelihood for the parameter uq (in units 



of 10" 



while using the polynomial of order 5 



pansion, we ended up with a marginal like lihood peaked 
aroun d qo « 0. Thus whereas the data in (jKnop et. al. 



2003) validated the claim of accelerated exp a nsion , the 
Constitution SN dataset in Ifficken et. ah] ||2009j) fa- 
vors a coasting evolution; i.e., the universe may neither 
be accelerating nor decelerating. However, the presence 
of substantial amount of dark energy and dark matter 
would still be required to explain the data. 

Here one observes that the result qo ~ could be 
connected to the value of h and also that properly 
including h in the analysis may further decrease the 
constraining power of the data. Therefore one should 
explore the consequences of using Gaussian priors on 
h from other measurements too. In fact, we have 
marginalised the likelihoods over the Hubble parame- 
ter, with a Gaussian prior h = 0.68±0.06, as mentioned 
above. But since the likelihood curve for h obtained 



F ig. 5 Margi nal likelihood for the parameter qo, obtained 
m I John I (l2005h. while us i ng the polynomial of order 5 and 



the data in 



Knop et. al.~l {2003). 



from the previous analysis by John ( 20051 ) is not very 
sensitive to its variation in the concerned range (unlike 
the case of go, r o, etc.), it is more appropriate to em- 
ploy priors for h deduced from other measurements. We 
propose that this procedure shall be followed in future 
analyses. 

The considerable spread left in the marginal like- 
lihoods shows that even now there is some freedom 
in choosing the values of those parameters. In other 
words, there is a sizable volume in the parameter space 
that can have the same low \ 2 ■ But this should not be 
viewed as a drawback of the analysis; instead, this sim- 
ply reflects the fact that the data are not yet accurate 
enough. Some recent analyses of Constitut ion SN data 



endorses this result (jShafieloo et.al. 112009). This free 



dom in SN data was noted earlier ( John 



2004 120051) 



which highlights the strength of the Bayesian model- 
independent approach. 

Based on the mean values obtained for these param- 
eters, we compute the successive terms in the series @. 
With time in units of 10 17 s, the series can be written 
as 



1 + 2.106 x 10 _1 T - 2.22 x 10~ 2 <? T 2 



+1.55 x 10 _3 r T 3 



0.819 x 10" 4 soT 4 



-3.45 x 10~ 6 u T 5 



where we have taken h — 0.65 (only to evaluate this 
series) . With the values of the parameter in the ranges 
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obtained in the analysis, this series appears to converge 
even for | T | as large as ~ 3 x 10 17 s. However, this 
feature is not essential for our analysis, for we have 
assumed only a polynomial form for the scale factor. 
The situation was not different in the previous work 
cither. 



3 Conclusion 

We assumed that a Taylor series form for the scale 
factor a(t) is valid and attempted to find the coef- 
ficients in this expansion using the recent Constitu- 
tion SN data. The new marginal likelihoods obtained 
for its coefficients give valuable information regard- 
ing the expansion history of the universe. It is found 
that there is significant constraining of these param- 
eters wh en compared to pre vious analyses using the 
data in ( Knop et. al. 20031) . The shift in the com- 



puted mean value of the deceleration parameter qo, 
from that found in the previous analysis is noteworthy. 
Even when we start with a prior probability distribu- 
tion that strongly favors an accelerating universe, the 
marginal likelihood for the deceleration parameter ob- 
tained from the present analysis using the Constitution 
data is found peaked around qg — 0. However, we re- 
iterate that the considerable spread still found in the 
likelihoods of these parameters indicate freedom in the 
choice of their numerical values. 

A distinguishing feature of our analysis is that the 
marginal likelihoods for each parameter obtained in the 
previous case is chosen as the prior probability distri- 
bution in the present one, thereby implementing the 
Bayesian method in true spirits. The work is also in- 
tended as a demonstration of this fundamental require- 
ment in Bayesian analysis. However, we have noted 
that the results obtained in this paper may heavily de- 
pend on the prior chosen for h. Thus it is important 
to evaluate expansion rates using prior for h deduced 
from other measurements too. It is expected that in 
future when the SN dataset becomes large enough, the 
expansion coefficients get sharply peaked marginal like- 
lihoods and become the most basic model-independent 
description of the expansion history of the universe. 
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